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ABSTRACT 

The early evolution of planetary systems is expected to depend on various periods of disk 
matter accretion onto the central star, which may include the accretion of metal-rich matter 
after the star settles on the main sequence. When this happens, the accreted material is rapidly 
mixed within the surface convective zone and induces an inverse mean-molecular- weight gradient, 
unstable for thermohaline convection. The induced mixing, which dilutes the metal excess, may 
also have important consequences on the light elements abundances. We model and analyze 
this process, and present the results according to various possible accretion scenarios. We give a 
detailed discussion of the different ways of treating thermohaline mixing, as proposed by previous 
authors, and we converge on a consistent view, including the most recent numerical simulations. 
Wc show how the observations of light elements in stars can be used as tracers of such events. 

Subject headings: starsiabundances, stars: solar-type, stars: convection, exoplanets, planetary systems 



1. Introduction 

Solar-type stars are currently divided into 
exoplanets-host stars (EHS), which are known 
to harbor planets, and non-exoplanets-host stars 
(NEHS), for which no planet has yet been detected 
with the available techniques and instruments. 
This dichotomy is obviously related to exoplanet 
detection bias, as stars without detected planets 
may very well host planets in the presently non- 
detectable range. In this framework, the studies 
concerning the observed differences between these 
two categories of solar-type stars must be taken 
with caution. However, taking this restriction 
into account, these studies give interesting infor- 
mation, at least on the observed EHS. 

Here we present a theoretical study concern- 
ing the fate of the metal-rich (planetary) material 
which may fall onto EHS during their early lives, 
and we analyze the consequences of such events 



on the abundances of light elements. We compare 
our results with the presently available observa- 
tions, and give some predictions and discussion for 
future work. 

The observed EH S are, on average, metal-rich 
compared to NEHS (jSantos et al.luOOll) , although 
some v ery metal poor EHS have also been ob- 
served ( Setiawan et al. I 20101 ). This average metal 
excess was a subject of discussion during the last 
decade, as the star could either has started with a 
higher metal content than average, or be related 
to metal-rich accretion onto the star. This last 
explanation is now ruled out for several reasons 
including the fact that the metal-rich accreted 
matter cannot remain inside the outer convective 
zone but falls down inside the star due to ther- 
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mohaline convection. IVauclaiij ()2004|) gave a first 
discussion of this specific process. She showed 
how the accreted matter was mixed down inside 
the star and left only a very small /x-gradient at 
the end of the mixing process. It is now proved 
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that this remaining ^-gradient is too small to 
account for the o b serve d overmetallicity. Con- 
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trary to IVauclaiil (|2004l ). 
suggested that the induced mixing region could 
go down t o the lithium destruction zone. Inde- 
pendently, [Caraud ( 2011 ) performed precise com- 
putations of thermohaline convection in a static 
IAMq model, with the assumption of accretion of 
one Jupiter plan et, using the thermo haline cocfB- 
cient derived by iTraxler et al. I (|201ll) (see below, 
section 5.7). She confirmed that the remaining 
metal excess was much too small to account for 
the observed overabundances, and found that the 
mixing region extended below the lithium burning 
region. 

The cfhcicncy of the thermohaline mixing and 
the resulting abundances strongly depend on the 
atomic diffusion processes at work in the stars. For 
example in the case of Carbon Enhanced Metal 
Poor Stars (CEMP), which accrete matter while 
they have already spent some time on the main 
sequence, the thermohaline mixing is strongly lim- 
ited by the pre-exis ting stable M-gradient, u iduced 
by helium settling ( Thompson et al. 20081 ). This 
is not expected to occur in the case of EHS which 
accrete matter at the very beginning of their lives, 
when atomic diffusion did not yet have time to 
built important helium gradients. However, as will 
be seen below, even in this case atomic diffusion 
is not negligible and has to be taken into account. 

Wc precisely analyze the consequences of metal- 
rich accretion onto solar-type stars, for various 
stellar masses and various assumptions of accreted 
masses. We simulate the induced thermohaline 
mixing using various prescriptions for the mixing 
coefficie nt, including t he mo st recent 2D simula- 
tio ns by Denissenkovl (2010l ) and 3D simulations 
by iTraxler et al.l ( 2011 ) and we compare the re- 
sults. We show that the situation may be quite 
complex. In any case, if some metal-rich layer can 
remain in the stellar outer regions, it is much too 
small to be observable. This confirms the primor- 
dial origin of the average overmetallicity observed 
in EHS. On the other hand, lithium and beryl- 
lium destruction layers are reached in most cases 
and these elements may be partially destroyed. 

In the present computations, we do not add 
any other mixing process which can also destroy 
lithium and beryllium, like rotation-induced mix- 
ing or mixing by internal waves. The reason is that 



wc want to disentangle the various processes act- 
ing on the light elements and look at the effects 
of thermohaline convection only. In the future, 
these other mixing processes should be added for 
a complete description of light element evolution 
in solar-type stars. 

In section 2, we give a review of the possible 
accretion phenomena expected on EHS and their 
observed evidences; we also discuss in this sec- 
tion the induced thermohaline process. A review 
of the observed element abundances in EHS ver- 
sus NEHS stars is given in section 3. Section 4 
is devoted to the modeling techniques. Section 5 
presents our computations and results; the effects 
of the various parameters like stellar mass, accre- 
tion mass, different prescriptions for the mixing, 
etc. arc compared in detail in this section. Finally 
the consequences of these results arc discussed in 
section 6, followed by the conclusion. 

2. The accretion phenomenon and the 
induced-instabilities 

2.1. Accretion onto exoplanet host stars 

2.1.1. The planetary system evolution 

It is now well-established that stars may ex- 
perience episodic accretion along their evolution 
through different processes. The planet formation 
period which extends from the star b irth to the 
end of the PMS phase (jlda fc Linll2008l ) is a favor- 
able epoch for such events. The interactions of a 
young planet with its surrounding disk affect the 
planet's orbital energy and angular momentum. 
As a consequence the planet may migrat e radially 
towa r ds or away from th e central star (jTerguem 
20101 : iLubow fc Idall2010l) . 

At the end of the protoplanetary disk phase, 
the planetary system architecture is not neces- 
sarily set. Current models for planetary system 
evolution indicate that various instabilities may 
commonly arise within the first gigayear after the 
system birth, as it moves towards a stable con- 
figuration. During this phase, the system experi- 
ences quasi-stable periods interrupted by more or 
less violent restructuring episodes. These events 
which may result in orbital reconfigurations, colli- 
sions, gravitational captures or ejections, may be 
favoured and intensified by the presence of aster- 
oids and debris belts. During these catastrophic 
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events, planetesimals, asteroids or debris may be 
pushed towards the central region and fall onto 
the inner objects and/or the central star. 

2.1.2. The Solar case evidences 

Evidences for the previously described evolu- 
tion are present in our Solar System and comes 
from a variety of sources. 

Planetary formation theories suggest that the 
giant planets formed on circular and copla- 
nar orbits different from their current path 
(jTsiganis et al.l 120051) . The asteroid and Kuiper 
belts are believed to have been formed with con- 
siderably mo r e mass than their present c onten t 
(|Sternl Il996l : lO'Brien. David and Bott"k3 l2007t) . 
their radial lo c ation may have changed over time 



(Gomes et al.l l2005t Thommes. Brvden fc Rasiol 



20071 ) and the dynamical structure of the Kuiper 
belt provid es evidence for scu l pting by planet 



migra tion (jChiang et al.l l2007t iMorbidelli et al 
20081 ). The petrology and cratering records on 
the Moon and on telluric planets suggest that a 
cataclysmic spike in the cratering rate occurred 
nearly 700 million year s after the planets formed 
(jHartmann et al.ll2000l) . This event which is re- 



ferred to as the Late Heavy Bombardment (LHB). 
is characterized by an intense period of planetes- 
imal bombardment pro bably originated fr o m the 



existing asteroid belt ( Kring et al. 2002; Tagld 



I2OO8I ) a nd spreading throu ghout the inner Solar 
System (jStrom et al.ll2005l ). This episode would 
have last 50 to 200 million years. 

Models and theories ha ve been proposed to 
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explain these observations (iTsieanis et al 
Thommes. Brvden k, Rasio 2007t Chiang et al 



the components of the asteroid belt would have 
been ejected from the planetary system or simply 
pushed towards outer regions (forming the current 
Kuiper belt). However some objects of the belt 
would have experiment high orbital eccentricity 
motions allowing them to enter the inner Solar 
System and impact with telluric objects, resulting 
in the LHB (|Gomes et al.ll2005t ). 

The solar system history probably reveals a typ- 
ical scenario for planetary system evolution. Var- 
ious models indicate that the instabilities encoun- 
tered by the Solar System may be a common oc- 
currence for planetary systems and that this can 
occur after up to 1 Gyr of quasi stable evolution. 

2.1.3. Interference by asteroid belts 

The presence of solid bodies belts may inter- 
fere in the stabilization process of planetary sys- 
tems. In particular, as stated in the previous 
section for the Solar System case, perturbations 
of asteroid and/or debris disks may result in the 
scattering of planetesimals throughout the plan- 
etary system, favouring the occurrence of intense 
impact episodes. As most stars are supposed to 
form surrounded by protoplanetary disks, they 
are expected to present remnant belts (seen as 
asteroids or debris disks) resulting from a partial 
clearing process of their original disks. Indeed 
massive circumstellar disks are observed around 
several hundred main sequence stars and preferen- 
tially around exoplanet host stars (hereafter EHS) 
( Brvden et al] l2009t IWvatt. Clarke and Grea^ 



20071 ). We consequently expect planetary systems 



20071: iFord fc Chiang"200f:'Thommcs'Dunca.n and 
Gomes et alEoOS : Morbidclh ct al. 2005: Tsiganis 
20051 ). According to the most popular ones, at the 
end of the protoplanetary disk phase, the gas giant 
planets were surrounded by a rich and interacting 
asteroid belt imposing slow planet migrations. As 
a consequence of these migrations, the occurrence 
of a 1:2 orbital resonance reached by Jupiter and 
Saturn would have been the catalyst for a catas- 
trophic episode. This phenomenon would have 
lead to a dramatic destabilization and a rapid 
restructuring of the outer Solar System, involv- 
ing the giant planets and the asteroid belt. As 
a result of this process, the giant planets would 
have evolved to their current orbits while most of 



to offer preferential conditions for violent bom- 
bardment episodes and consequently for accretion 



Levison event onto the central star. 
et al.l 

2.2. The thermohaline instability 

The accretion of planetary (metal rich) material 
onto a star produces a brutal increase of its sur- 
face metal content. The metal-rich accreted ma- 
terial is rapidly mixed within the surface convec- 
tive zone: a metal-enriched zone may then overlay 
lighter layers which produces an unstable weight 
gradient in the external layers of the star. A simi- 
lar situation is often encountered in oceanography, 
when warm salted water lies upon cool unsalted 
one. Large inverse salinity gradients make warm 
salted blobs diffuse rapidly downward in spite of 
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the stabilizing temperature gradient. Heat diffus- 
ing more rapidly than salt, the warm salted blob 
falling down in cool fresh water experiences a tem- 
perature decrease before the salt has time to dif- 
fuse out; the blob continues falling because of its 
weight, while unsalted matter rises around. This 
special mixing process which takes the form of 
elongated cells, the "salt fingers", is called "ther- 
mohaline convection" . It proceeds until the fingers 
mix with the surroundings and a stable salinity 
gradient is restored. 

A similar type of convection can grow inside 
stars. It occurs in a stable stratification that sat- 
isfies the Ledoux criterion for convective stability: 



Vad - V™d + > 



(1) 



but where the molecular weight decreases with 
depth: 



(2) 



With the classical notations, Vad and V are re- 
spectively the usual adiabatic and local gradi- 
ents (dlnT/dlnP), = {d\n p/d\n ^)p^t and 
S=id\np/d\nT)p^f,. 

In the stellar case, the role of the aquatic salin- 
ity is then played by the molecular weight. When 
a blob of heavy matter (i.e with high molecular 
weight) begins to fall down, heat is exchanged 
with the surroundings more rapidly than parti- 
cles so that it goes on falling like a finger, forc- 
ing light matter to go up instead. This hap- 
pens every time heavy element accumulation oc- 
curs above lighter layers. It leads to short episodes 
of deep thermohaline mixing which proceed until 
the molecular-weight gradient becomes too small 
to trigger the instability. This ca n be summarized 
wit h the following condition (see IVauclaiii (j2004l ) 
and Traxler et al 



Lg conditL 

Jlmil)): 



1 < i?o < 



1 



where Rq represents the ratio: 

j-> ^ ad ^ rad 

^0 = 



(3) 



(4) 



and r is the inverse Lewis number, ratio of the 
molecular diffusivity to the thermal diffusivity. 



A more complete description of the methods 
used to model thermohaline convection is given in 
Section 4.2. 

Thermohaline instabilities due to inverse fi- 
gradients are expected to arise in stellar interiors 
in a variety of already identified contexts. This is 
the case for instance when helium or carbon rich 
material is deposited at t he surface of a star in a 
mass-transferriiig binary ( Stothers fc Simonlll969l : 
Stancliffe et"all l2007t [Thompson et al.l |2008[) . It 
also develops in stars where radiative levita- 
tion pr oduces heavy element accumulation re- 



gions (|Theado et all l2009l ). The thermoha- 



line mixing has also been proposed as a possi- 
ble mechanism to account for the photospheric 
com position of low and intermediate mass gi- 
ant ('Eeaiet on. Dearborn fc La ttanzio"2006! . l2008t 
Char b onnel fc ZahnI l2007t ICanticllo fc Langer 
20081: ISiessI 120091 : ICharbonnel fc Lagardd 120101: 
Stancliffd l2010l) . In such stars an inverse /i- 
gradient is created by the ^He('^He,2p)^He reac- 
tion in the external wing of the hydrogen burning 
shell. An inverse /i-gradient may also arise after 
the ignition of ^He burni ng in a degenerate sh ell 
(i.e. core helium flash; see iThomasI (Il967lll970l) '). 

Up to now the effects of accretion-induced ther- 
mohaline instabilities on EHS have b een little 
studie d. They have been di scussed bv | Vauclair 
(|2004h and more recently by iGaraudI (|201ll ). but 
no detailed stellar evolutionary models have yet 
been computed to investigate the effects of accre- 
tion and thermohaline convection on the stellar 
lithium, which is the aim of the present paper. 

3. Observational background : abun- 
dances in exoplanet host stars 

3.1. The overmetallicity of EHS 

Stars hosting planets compared with other stars 
of the same spectral types clearly show a metal- 
licity excess of a factor 2 on average. Two sce- 
narios have been proposed to explain the metal 
enrichment of EHS. The first one, the accretion 
hypothesis, considers that EHS formed with nor- 
mal abundances but experienced during the planet 
formation period metal-rich mater ial accretion 
(|Gonzalej[l998t iMurrav et al.ll200l[ ). The inges- 
tion of planets, planetesimals, asteroids, comets, 
dust would then produce a local surface over- 
metallicity. The second scenario, the primordial 
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hypothesis, assumes that central stars of plan- 
etary syste m form out from metal rich molecu- 



(|2006l ) but cont radict by iLuck fc Heiteil (l2006l) 



lar clouds ( Pinsonneault et alJl2001: Santos et al 



200 ll 120031 ) : the overmetallicity of the protostel- 



lar gas would be a necessary condition for planets 
to form around stars. In this case, the stars would 
be metal-rich from the center to their surface. 

Several studies have tried to differentiate be- 
tween the two possible explanations using dif- 
ferent techniques i nclud ing asteroseismology (e.g. 



Bazot fc Vauclairl 120041) . Today a consensus 



seems to be achieved which favors the primor- 
dial hypothesis. Several arguments have been 
invoked which rule out the accretion scenario. 
It was first thought that the dilution of metal- 
rich material in stellar external convcctivc zones 
should lead to metallicity variations depending 
on the stellar mas s , whi c h were not observe d 
(IPinsonneault et al.l l200lt ISantos et ahl l200l[ ). 



Another argument was that the overmetallicity 
observed in EHS would require the ingestion of 
several tens of earth, which is not completely ex - 
cluded but seems quite high (jlsraelian et al ]|2004) . 
More recently, as already discussed in previous sec- 
tion, it was recognized that the ovcrmctallic mat- 
ter diluted in the surface convcctivc zone would 
rapidly sink in the radiative interior due to thcr- 
mohaline instability effects, so that only a small 
fraction of the accr eted matter cou ld re main in 
the ou ter layers (see IVauclaii (2004) and GaraudI 

(tmH)). 



3.2. The lithium and beryllium enigma 

3.2.1. Observational context 

The potential discrepancy between lithium 
abundances in EHS and NEHS is a controversial 
subject which has divided the observers commu- 

Lawd 



Gonzalez 



nity f or more than 10 years. 

()2000l) irst suggested that EHS present smaller 
hthium abundances than fi e ld sta rs. They were 
soon contradic t by RvaiJ ( 2000f ) and later by 
Gonzalez et al. ( 200ll ) who stated that lithium 
abundances of stars with detected exoplanets are 
indistinguishable from those of s tars w ithout de- 
tected planets. Ilsraelian et al. (l2004l) . who re- 



visited this topic, reported a lithium deficiency 
in EHS within the effective temperature range 
5600K-5850K. Their results were confirmed by 
Takeda fc Kawanomotol ( 2005 ) and Chen fc Zhaol 



More recently 
(j2008( ) and Ilsrael ian et al 



Gonzalea 



Takeda ct al.' ('2007|), 

(2009) confirmed an 
enhanced lithium depletion in Sun like stars 
with orbiting planets bu t once more th e se re - 
sults were contradict by Baumann et al. ( 2010l ) 



who attributed the l i thiura depletion derived by 
Israelian et aD (|2004l l2009l ) to a wrong interpre- 
tation of observational data or to observational 
bias. 

The situation is sti ll more diffic ult for beryl 



hum. Delgado-Menal ( 201l[ ) and Galvez-Ortiz 



(|201ll ) have determined the beryllium abundances 
in a sample of solar-type stars and tried to de- 
termine a difference between cxoplanct-host stars 
and stars without detected planets. They may find 
an extra beryllium depletion for the coolest stars 
(Te//< 5500 K), but this result stih has to be 
confirmed. 

3.2.2. Theoretical context 

Considering the potential lithium deficiency ob- 
served in EHS as true, theoreticians proposed sev- 
eral explanations to account for this "extra" de- 
pletion. Some of them suggest that the lithium 
destruction and the planet formation both results 
from special conditions surrounding some stars 
(e.g. overmetallicity, long-lived disk) . Others con- 
sider that the destruction is directly related to the 
presence of planets. 



Metallicity effects. ICastro et all (|2009l ) 



puted stellar evolutionary models including atomic 
diffusion and rotation-induced mixing. In this 
framework they investigated the effects of var- 
ious primordial overmetallicities on the lithium 
depletion. They showed that the EHS-lithium 
deficiency cannot be attributed to overmetallicity 
effects : extra-mixing is required below the surface 
convective zone to account for the Li observations. 



Stellar rotational history. iBouvierl (|2008l ) in- 
vestigated a possible connection between a lithium 
deficiency, the presence of planets and the rota- 
tional history of solar type stars. Thanks to rota- 
tional evolution models, he showed that contrary 
to ZAMS fast rotating stars, the slow ZAMS ro- 
tators develop strong differential rotation between 
their surface convective zone and their radiative 
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interior, which leads to hydrodynamical instabil- 
ities responsible for enhanced lithium depletion. 
He concluded that EHS were slow rotators on the 
ZAMS and that the protoplanetary long-lived disk 
probably responsible for their slow rotation rate 
could also be a necessary condition for the planet 
formation. This could explain the correlation ob- 
served between the presence of planets and large 
lithium destructions. 

Formation, migration and accretion of plan- 
ets. Several authors suggest that the strong 
lithium depletion observed in EHS could be di- 
rectly related to the mechanism of formation or 
migration of pla nets. This idea is sup ported by ob- 
servations from Chen fc Zhaol ( 2006[ ) which show 
that overmetallic stars without detected planet do 
not present strong lithium depletion. 

Israelian et al.l ( 20041 ) suggest that protoplane- 



tary disks lock a large amount of angular momen- 
tum during the planet formation mechanism thus 
inducing some rotational braking in the host star 
during the PMS. The mixing resulting from the 
induced differenti al rotation could lead to a large 
Uthium depletion. Ilsraelian et al.l (j2004l ) also sug- 
gests that the lithium destruction could be asso- 
ciated with a late migration of giant planets at 
the end of the PMS. Angular momentum trans- 
fer from the planetary disk to the outer layers of 
the star and/or tidal effect of the giant planet on 
the host star could induce differential rotation and 
then mixing in the host stars. 

If the migration process is overly; efficient, plan- 

ets ma y be accreted onto the star. iMontalba n fc Rebolo 



(|2002l ) investigate the effects of the accretion of 
Jupiter-like planets on the lithium isotopic abun- 
dances in main sequence stars, using either stan- 
dard stellar models or models with a light mixing 
process induced by gravity waves. They do not 
take into account any thermohaline convection, so 
that the stellar outer layers are not mixed with the 
lithium depleted regions. As a consequence, they 
obviously find an increase of the lithium abun- 
dance and the 6/7 isotopic ratio due to the planet 
chemical composition. Baraffe & Chabrier (2010) 
explores the effects of early episodic accretion (i.e. 
occurring during the PMS) on the internal struc- 
ture of young (fully or nearly fully convcctive) low 
mass and solar type stars. They conclude that 
an early protostar accretion history can affect the 



internal structure of the star and can severely 
enhanced the Li-depletion in these object. 

4. Modeling 

We computed stellar models of solar type 
stars, using the Toulouse- Geneva Evolution Code 
(TGEC). T he code is described in detail in 



Richard et al. 



(|2004) and 
All our models include atomic 
desc ribed 
()2009l) 

number of individual elements 



puted 



as 



Theado et al 



Hui-Bon-Hoal (|2008f ). 

diff usion com- 

Hui-Bon-Hoal (|2008l ) and 
The selective diffusion of 



a large number ot individual elements impor- 
tant for the stellar structure is computed with 
respect to hydrogen, including the effects of 
the local pressure gradient, temperature gradi- 
ent, concentration gradient and radiative trans- 
fer. At each time step, all the abundances are 
renormalized to account for the correlated mo- 
tion of hydrogen and the induced modifications 
of the opacities are taken into account. When 
macroscopic motions occur, which can be de- 
scribed in terms of a macroscopic diffusion co- 
efficient, the mixing effect is added as a new 
term in the diffusion equation, acting on the 



^ ^ ^^^^^v.^^.., v^^ ^^^^^^ ^ 

pers (e.g. Vauclair. Vauclair and Michaud 




1978h, 


IVauclair & Vauclaii-' (■ 19821). Richard et al. 




19961). 


iVauclair & Theado (2003i). Michaud. Richer & RichardI 



In the TGEC code the main isotopes of the fol- 
lowing elements are treated separately: He, Li, Be, 
B, C, N, O, Ne, Mg, Ca, Fe. The other species are 
treated together. For each separate isotope i, the 
diffusion equation is written; 



d{pci) 
dt 



div{pCiVi) = 



(5) 



where c; is the concentration of the isotope and 
(pciVi) the diffusion flux. In atomic diffusion the- 
ories, Vi is traditionally called the "diffusion ve- 
locity". The expression is deduced from detailed 
theories of transport proc esses in multicompo- 
nent non-uniform gases (e.g. Chapman fc Cowlind 
1970f ). For test atoms, it may be written: 



Vi = D, { _£l±_El!Ly In c + fcpV In P 



D, 



-FfcrVlnT- 



kT 



(6) 
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where the first term is the classical diffusion term 
and the three other ones are respectively the pres- 
sure, ther mal and radiative diffusion terms (e.g. 
Eq. 2 of iTheado fc Vauclai j [2003al) . D, is the 
atomic diffusion coefficient and Dth the mixing dif- 
fusion coefficient. In the cases where atomic dif- 
fusion is negligible, these equations reduce to the 
classical diffusion (mixing) equation. 

In the following, we discuss how accretion and 
thermohaline convection are introduced in the 
code. 

4.1. Accretion 

We assume that stars could swallow planetary 
material during their early main sequence life. The 
accretion of a given amount of matter is modeled 
through an increase of the metal content of the 
surface convective region. The accreted mass is ex- 
pressed in Jupiter Mass (hereafter Mjup) or Earth 
Mass unit (hereafter M^); its composition is as- 
sumed similar to the sola r mixture determined by 
Grevesse fc Noelsl (|l993h but devoid of hydrogen 
and helium. The composition of the surface con- 
vective zone is computed again after each accretion 
episode : the mass fractions of the chemical ele- 
ments are computed taking into account the newly 
accreted matter but the structure of the star is 
supposed unaltered by the accreted material. 

In the case of large amounts of accreted matter, 
the /i-gradient is large enough at the beginning 
to destabilize the medium by dynamical convec- 
tion. This is modelled by a very large diffusion 
coefficient simulating complete mixing. Thermo- 
haline mixing begins when the /i-gradient is small 
enough for the medium to be stable against dy- 
namical convection. The prescriptions used to de- 
scribe this mixing are discussed below. 

4.2. Thermohaline convection 

The effects of thermohaline convection as a mix- 



Ulrich 



i ng pr oces s in stars are far from t rivial 

(|l972l ) and lKippenhahn et al.l (|l980l ) (hereafter re- 
ferred to as KRT) proposed to model the thermo- 
haline mixing process in terms of a diffusion coef 
ficient proportional to the inverse /i- gradient: 



AacT^ Hp I d In /i I 
icpnp^ Vad - V dr 



(7) 



which can simply be written: 

Dth = CtKrUa ^ 



(8) 



where kt is thermal diffusivity, Hp is the pres- 
sure scale height, a the radiation constant, c the 
speed of light, T the temperature, Wad the adi- 
abatic gradient, cp the specific heat at constant 
pressure, n the opacity. The variable Ct is a fac- 
tor involving the aspec t ratio (lengt h/ width) of 
the fingers. KRT and lUlrichI |l972[ ) considered 
different geometries for the fingers and, as a con- 
sequence, they de duced d if ferent Ct values (12 for 
KRT and 658 for lUhichI (Il972l)). In their com- 
putati ons of red giant stars, ICharbonnel fc Zahn 
()2007l) considered this parameter as free and ad- 
justed it {Ct = 1000) to reproduce the abundance 
observations in evol ved stars. On the other hand 
iTheado et aD (|2009[ ) used the original KRT value 
with Ct = 12 for the computations of diffusion 
induced thermohaline convection. 

A large step forward has recently been acom- 
plished o wing to the numer ical simul ations per- 
forme d bv lPenissenkov ( 2010[ ) (2D) and lTraxler et al 
(|201l[ ) (3D). They both deduced from their simu- 
lations of fingering convection a near-unit aspect 
ratio of the fingers much sm aller than that used by 
Charbonnel fc ZahnI (2007). In the intermediate 



regime, their different expressions for the thermo- 
haline mixing coefficients lead to values very close 
to the KRT ones, if computed with Ct = 12. Im- 
portant differences appear only for the boundary 
values of the Rq ratio, when Rq — > l/r, which cor- 
responds to the case of very small /^-gradients, or 
i?o — > 1, which corresponds to the largest possible 
/i-gradients before dynamical convection. 

The Denissenkovl ( 2010h coefficient may be writ- 
ten as: 



Dth = CtKTiRo - ly'i^ - Rot) 



(9) 



which reduces to the KRT one for 1 << i?o << 
l/r, but vanishes for Rq = 1/t and goes to infinity 
when i?o = 1- 

Traxler et all (|201l[) (hereafter TGS) have per- 



formed 3D simulations of thermohaline convec- 
tion with various sets of physical parameters, and 
found that the ratio of the turbulent to the diffu- 
sive fluxes settled onto a single universal profllc for 
all parameter values, which seems a robust result. 
They derive a completely empirical expression for 
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the mixing coefficient which, at first sight, looks 
quite different than the previous ones, but present 
also some similarities: 



Dth = lOly^K^exp (-3.6r)(l - r) 



1.1 



(10) 



where, in stellar conditions for which r ^ 1, 
r = {Rq — 1)t. After some manipulations, we find 
that in the intermediate regime this mixing co- 
efficient may be approximated as Dth = Ctktt, 
which is slightly smaller than the KRT one with Ct 
= 12. Once again, important differences only ap- 
pear for the behavior in the extreme cases. Here 
Dth vanishes for Rq = 1/r but it reaches a fi- 
nite value when Rq ~ 1. This last point will be 
discussed in section 5.7. An illustration of this 
discussion may be found in [14] where we present 
a graph of the KRT and TGS mixing coefficients 
precisely computed in one of our models. 

In the following, we first present results ob- 
tained with the KRT prescription using Ct = 12, 
and we also do computations with the assump- 
tion Ct ~ 1000 for comparison purposes. We 
show that in this case, the resulting lithium and 
beryllium depletions are too high to be realistic, 
which confirms that this coefficient is overesti- 
mated. Then we present computations obtained 
with the new TGS empirical expression and we 
compare it with the previously obtained results. 

5. Computations 

We made computations for stars with masses 
between 0.70 and 1.30Mq. We tested various ini- 
tial metal contents with metallicity values rang- 
ing from [Fe/H]=-0.20 to +0.20. We computed 
models with a solar initial helium mass fraction 
and models with an initial helium abundance vary- 
ing with metallicity (Yo=Yg) according to the 
galax y chemical evolution law (jlzotov fc Thuan 
One to fifteen accretion events were intro- 
duced in our models with various accretion rates 
for a total accreted mass ranging from IMg to 
l.SMjup. We present here a choice of our models, 
representative for all the computations done. 

5.1. The accretion/thermohaline mixing 
sequence 

In this section we describe in detail the accre- 
tion/thermohaline mixing process, and present re- 
sults obtained with the KRT coefficient {Ct = 12). 



5.1.1. The case oj one accretion episode 

As a first example, we present the result ob- 
tained for a 1.10 M0-model with Yo=Yg=0.29 
and [Fe/H]=0.20. This model experiences a sin- 
gle accretion event of O.OSMjup two million years 
after its arrival on the ZAMS. This accretion rate 
is of the order of a giant planet core, which repre- 
sents a realistic value for the mass of heavy mate- 
rial collected by the star through an impact with 
a Jupiter like object. 

Figure [T^ displays the molecular weight pro- 
file in the outer regions of the model at vari- 
ous epochs during the "accretion/mixing period" 
(which refers to the period starting with the ac- 
cretion episode and ending when the thermohaline 
mixing stops). 

The instantaneous dilution of a O.OSMjup plan- 
etesimal in the external convective zone produces 
an enhancement of the surface metal content of 
the model. The dotted (purple) line on Figure [1^ 
represents the theoretical metallicity increase pro- 
duced by the dilution in the convective zone of the 
impacting object. This metallicity increase cre- 
ates a sharp positive /i-gradient at the transition 
between the radiative interior and the convective 
envelope which triggers thermohaline instabilities. 
The induced mixing acts to soften the unstable /i- 
stratification. It occurs simultaneously with the 
atomic diffusion and in particular the downward 
migration of helium which also contributes to re- 
store a stable /i-gradient. When a flat /i-gradient 
is reached below the convective zone, thermoha- 
line convection becomes inefficient and a stable /i- 
gradient induced by helium settling begins to take 
over. Thermohaline convection still occurs below 
this stable region. As will be seen in forthcom- 
ing sections, this effect is less important for larger 
accretion rates (section 5.3). 

Figure [TJd displays the diffusion coefficient used 
to parametrized the thermohaline mixing at the 
same epochs as presented on Figure [1^. The mix- 
ing starts right after the accretion event, when 
an inverse /i-gradient is created at the transition 
between the radiative and convective zones. It 
mixes a localized region restricted to the unsta- 
ble /i-stratified layers lying below the convective 
zone. This process smooths down the unstable 
//-gradient in the mixed zone, it decreases the dif- 
fusion coefficient but it also continuously widens 
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the unstable region. The bottom of the mixed 
zone deepens until a stable stratification can be 
reached : the mixing then stops and the diffu- 
sion coefficient drops to zero. The phases during 
which the accretion mixing proceeds have strong 
consequences for light elements like Li and Be as 
discussed in the following sections. 

5.2. The multiple accretion event case 

Wc now examine the case of multiple succes- 
sive accretion events. We present the results ob- 
tained for two I.IOM0 (with the same Yq and 
[Fe/H] values as previously). In both cases the 
models experience ten accretion events. The first 
accretion episode occurs 2 Myr after the arrival on 
the ZAMS, and the following ones occur succes- 
sively every 2 Myr. In the first model, the amount 
of accreted matter is IMq (~0.003Mjup) at each 
episode. In the second model, the amount swal- 
lowed each time is O.OSMjup, of the order of a giant 
planet core, which represents a realistic value for 
the mass of heavy elements collected by the star 
through an impact with a Jupiter like object. 

Figure [H (a and b) and Figure |3] (a and b) 
present fj, and Dt profiles inside our models at 
various steps along the accretion/mixing period. 
As already discussed, the first accretion event in- 
duces a surface fi increase which triggers the thcr- 
mohaline instabilities. The resulting mixing tends 
to smooth down the unstable stratification before 
the next accretion event occurs. The next succes- 
sive accretion episodes continuously increase the 
inverse //-gradient, which feeds the mixing and in- 
creases the width of the mixed zone. After 2 ac- 
cretion events, the helium settling becomes unable 
to counteract and freeze the thermohaline mixing 
before the next event occurs. After the last accre- 
tion episode, the thermohaline convection still pro- 
ceeds, until a stable /x-gradient is restored. Dur- 
ing that process, a small stable region develops 
just below the convective zone. This is the reason 
why, on the top right graph, one curve, drawn af- 
ter 51 Myr, present a Dt which vanishes slightly 
below the convective zone. The time required to 
restore a stable stratification depends on the ac- 
cretion rate : the higher the accretion rate, the 
longer the mixing period and the larger the mixed 
region. The model experiencing the IMq accre- 
tion rate recovers a stable /^-stratification at 39 
Myr (i.e. 37 Myr after the first impact) while the 



model with the 0.03Mjup accretion rate is stable 
at 51 Myr (49 Myr after accretion beginning). 

As shown on Figures [2] and [3l the top of the 
Li-depleted region is close to the bottom of the 
surface convective zone. The thermohaline mix- 
ing then allows connecting the stellar surface with 
the lithium nuclear burning region and is there- 
fore expected to produce a Li depiction. Figures 
[5}: and[3J; represent the lithium profiles inside the 
two considered models. In both models the first 
accretion episode does not deplete significantly the 
lithium surface abundance: the mixing-induced 
lithium destruction is compensated or even ex- 
ceeded by the lithium amount brought by the im- 
pacting object. The next events increase the un- 
stable /i-gradient, as a result the diffusion coeffi- 
cient is enhanced and the mixed region extends 
towards deeper regions : as shown on Figures [2j: 
and [St, this leads to large lithium destructions. 
We will see in the next section that similar features 
are obtained for Be (the Be depletion is however 
expected to be smaller than the Li one since its 
nuclear burning region lies deeper in the stellar 
interior) . 

5.3. The accretion rate effects 

In this section, we test the effects of various 
accretion rates on the surface metallicity and the 
Li and Be depletions. In this framework. Figure 
m displays the variations of the surface metallic- 
ity over the accretion/mixing period for the two 
series of models (including the two models pre- 
sented in the previous section). All these models 
have the same mass (I.IOMq) and initial parame- 
ters (in particular the same chemical composition) 
but they experience various accretion scenarios : 
1, 5 or 10 accretion events of IM^ (0.003Mjup) for 
the first series of models and 1, 5 or 10 accretion 
events of 0.03Mjiip for the second series. 

Each accretion episode induces a metallicity in- 
crease which obviously depends on the accretion 
rate. For the considered accretion rates and num- 
ber of impacts, the overmetallicity does not exceed 
0.05 dex. After the end of the accretion period, 
the overmetallicity decreases under the combined 
effects of thermohaline mixing and atomic diffu- 
sion which makes heavy elements diffuse out of the 
surface convective zone. When the mixing stops 
(i.e. before 60 Myr for all the presented models), 
the remaining metallicity increase is smaller than 
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Fig. 1. — Description of one accretion/mixing period in a 1.10 model undergoing 0.03Mj„p accretion 
event 2 Myr after its arrival on the ZAMS. Figure a : molecular weight profiles in the external regions of 
the model. The dashed (black) line displays the atomic diffusion-induced /i-profile at 2 Myr just before the 
accretion episode, the dotted (purple) line represents the theoretical /i-profile resulting from the instanta- 
neous dilution of the planetesimal in the surface convective zone. The dashed-dotted (orange), long-dashed 
(magenta), solid (green) and long-dashed dashed (blue) lines represent the /t-profiles 2, 25, 54 and 89 Myr 
after the accretion event. After 4 Myr, the /^-gradient is still unstable below the convetive zone, whereas 
after 25 Myr, a diffusion-induced helium gradient leads to a stable /i-gradient in these layers, as discussed in 
the text. Figure b : profiles of the diffusion coefficient (cm^.s~^) modeling the thermohaline mixing (same 
line convention as in Figure a). The vertical line locate approximately the top of the Li- nuclear depletion 
zone. The grey regions locate the surface convective zone. After the accretion episode, mixing occurs deeper 
and deeper and links the outer layers down to the lithium burning region (dashed-dotted, orange line), until 
the competition with atomic diffusion begins to lead to stable layers below the convective zone (solid, green 
fine). After 89 Myr (long-dashed dashed, blue line), the thermohaline diffusion coefficient is flat and there 
is no more mixing. 
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Fig. 2. — a) molecular weight, b) thermohaline diffusion coefficient (cm^.s~^) and c) Li abundance profiles 
at various steps during the accretion/ mixing process in a I.IOMq model experiencing 10 accretion events of 
IMq. (Li/Lio) represents the ratio between the current lithium abundance and its initial value (i.e. at its 
arrival on the ZAMS, before accretion and diffusion start). The black lines represent the /i and Li-abundance 
profiles at 2Myr on the ZAMSjust before the first accretion event. The dark blue to magenta solid lines 
present the profiles at 4, 6, 8, 10 ... and 22 Myr. The light magenta dashed lines show profiles at 51 and 
97 Myr (the accretion has stopped but the thermohaline mixing still proceeds). For the two first accretion 
events, the same process as described in Figure 1 occurs. As an example, the horizontal arrow in the upper 
right panel locates the mixed region after 10 accretion events. The down right graph show the importance of 
this process on the surface lithium abundances. 
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Fig. 3. — Same figure as Figure [2] for a I.IOM© model experiencing 10 accretion events of O.OSMjup. The 
line convention is the same as in Figure [2] except for the light magenta dashed lines which show profiles at 
74, 167 and 236 Myr on the ZAMS. Here we can see that the small stable helium zone develops after 167 
Myr. 
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0.02 dex and is unable to explain the average over- 
mctallicity of 0.2 dex observed in EHS. 

Figure [5] shows the evolution of the Li and Be 
surface abundances over the accretion/mixing pe- 
riod. The results demonstrate that the accre- 
tion of planetesimals may decrease the surface Li 
and Be abundances, with a depletion rate depend- 
ing on the accretion rate and the number of im- 
pact events (the Li and Be depletion increasing 
with the accretion rate and/or the number of im- 
pacts). The Be nuclear depleted zone being deeper 
that the Li one, the obtained Be destructions are 
smaller than the Li ones. 

5.4. Influence of the stellar mass and the 
initial metal content 

In this section we present models with various 
masses and initial metallicities. All the presented 
models undergo 5 accretion episodes of O.OSMjup 
following the same scheme as previously discussed 
(same time intervalle between the ZAMS and the 
first impact and then between successive accretion 
episodes). As in previous sections, the thermoha- 
line mixing is still computed using the KRT coef- 
ficient with Ct ~ 12. 

Figure [6] presents the surface Li abundance dur- 
ing the accretion/mixing period in two sets of 
models computed with different initial metal con- 
tents. For each metallicity several masses are pre- 
sented. For a given metallicity, the Li-deplction 
induced by the combined effects of accretion and 
mixing decreases with increasing mass. Several 
competing effects lead to this result. 

The mass of the surface convective zone de- 
creases as the stellar mass increases, as a result the 
dilution of a metal rich planctesimal in the exter- 
nal convective zone should lead to larger unstable 
//-gradients and then to more efficient mixing and 
larger Li-depletion in more massive stars. However 
other competing effects lead to a different result. 
The Li amount brought through the ingestion of 
a metal-rich planetesimal increases the Li surface 
abundance of the target star. This increase is more 
significant in more massive stars (due to their nar- 
row convective zone). As a result more massive 
stars enter the mixing period with higher Li sur- 
face abundance. Moreover as the mass increases 
the top of the Li-dcpleted region moves away from 
the bottom of the convective zone. Larger mixed 
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Fig. 4. — Variations of surface ovcrmetallic- 
ity along the accretion/mixing period for several 
l.lOMg models. A[Fe/F] is the difference be- 
tween the current surface metallicity and its ini- 
tial value. Dashed lines are models experiencing 
(from bottom to top) : 1, 5 or 10 accretion events 
of IMg. Solid lines are models experiencing (from 
bottom to top) : 1, 5 or 10 accretion events of 
O.OSMjup. 



13 



3.5 - 



g.5 
< 



1.5 - 



1 



( 1 1 

- 


1 1 1 1 1 
1 


III! 

5xlM^ 


1 1 ) 1 1 1 1 1 1 1 

1 1 

lOxlMg ~ 


- V 






lx0.03Mj„p - 








5x0.03Mj„p : 








10x0.03Mj„p I 


1 1 1 


1 1 1 1 1 


1 1 1 1 


1 1 1 1 1 1 1 1 1 



1.6 



0.8 



I I I I 1 I I I I I I I I I I I I I ~ 



lxlM„ 5xlM„ lOxlM. 1X0.03M 




Jup 



5x0.03Mj„^ H 



10X0.03M 



Jup 



J I I I I I I I I I I I I I I I I I I I I I L 



_L 



50 100 150 200 250 50 100 150 200 250 
t (Myrs) t (Myrs) 

Fig. 5. — Variations of Li and Be surface abundances alon g the accretion/mixing pe riod for several I.IOMq 
models. The initial Li and Be abundances are taken from Grevesse fc Noel3 (1993) and respectively equal 
to 3.31 and 1.42. The curves are labelled according to the number of accretion episodes and the accretion 
rate undergone. 
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Fig. 6. — Lithium surface abundance over the accretion/mixing period in models experiencing 5 accretion 
episodes of 0.03Mjup. The presented models have different masses (0.8, 0.9, 1.0, 1.1, 1.2 and 1.3Mq) and 
initial metallicities. These models have been computed with the KRT coefficient, Ct~l2. 
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regions are then needed to connect the surface and 
the Li-nuclear burning zone. Both these effects 
(the accretion-induced Li increase and the grow- 
ing distance between the convective zone and the 
Li burning region) make the Li depletion smaller 
in more massive stars. 

Figures [7^ and b compare the Li depletion in 
models with the same mass but with different 
metallicities. It shows that the accretion/mixing 
induced Li depletion increases with the initial 
metal content. For a better understanding of this 
result we present on Figure [8] the surface con- 
vective mass as a function of the stellar mass 
for ZAMS models with various metallicities. It 
shows that the mass of the surface convective 
zone increases with metallicity. This feature con- 
firms that several competing effects leads to larger 
depletion in stars with larger surface convective 
zones. 



|9^). A retroactive effect of smaller /^-gradients is 
to decrease the diffusion coefficient (cf. Eq. [7]). 
The mixing efficiency then results from the com- 
petition between two opposite effects. The larger 
Cf tends to increase the diffusion coefficient and 
the mixing efficientcy, thereby leading to smaller 
/i- gradient, but on the other hand these smaller 
unstable gradients tend to decrease it. 

Figure [TU] illustrates the variations of the 
accretion-induced overmetallicity during the ac- 
cretion/mixing perio, using the large diffusion co- 
efficient (for comparison with Figure 4). Figure 
II II shows the lithium and beryllium depletion due 
to accretion and mixing. As previously discussed 
they depend on the accretion rate and the number 
of accretion events but for a given accretion sce- 
nario, the obtained destructions arc larger than 
those obtained with a smaller diffusion coefficient 
(cf. Figure [5]). 



5.5. Results obtained with the KRT coef- 
ficient and C't ^ 1000 

Up to now, we presented results obtained us- 
ing the prescription proposed by KRT with the 
adjustment Ct = 12. As already discussed, 
this coefficient is in much better agreement with 
the results of rece nt numerical computations 
(iTraxler et al. 2011 ) than the one proposed by 
Charbonnel fc ZahnI (l2007l) . using Q = 1000. We 
nevertheless computed models using this large 
value of Ct and give the results below for com- 
parison. 

Figure [9] present the results obtained with Ct = 
1000 for a I.IOMq model (with Yo=0.29 and 
[Fe/II]=0.20 as previously presented models), un- 
dergoing ten accretion episodes of O.OSMjup. The 
accretion sequence is the same as previously : it 
starts 2 Myr after the ZAMS beginning and the 
successive accretion events arc separated by 2 Myr. 

The models presented on Figures [3] and |9] differ 
only through the value used for the Ct factor. A 
larger Ct value increases the diffusion coefficient 
and consequently produces a larger mixing. The 
2 Myr-mixing occurring in-between two successive 
accretion events is more efficient in smoothing the 
unstable stratification. Consequently the succes- 
sion of several accretion episodes produces smaller 
unstable gradients in models with larger Ct values. 
This feature is indeed observed (cf. Figures|3K and 



5.6. Influence of the accretion events fre- 
quency 

In previous sections, we discussed either the 
cases of single accretion events, or the cases of sev- 
eral accretion events occuring every 2 Myr. When 
a new event occurs, the element abundances are 
first modified inside the convective zone due to 
the chemical composition of the accreted matter, 
and then they go on varying due to the thermo- 
haline mixing. In the case of lithium, each new 
event leads to a rapid lithium increase followed by 
a slower decrease induced by thermohaline con- 
vection and nuclear burning. Obviously the final 
results depend on the amount of accreted matter 
but also on the time scale between two events. If 
another accretion occurs before the complete set- 
tling of the previous one, the effects overlap. 

In this section, we discuss the effects on lithium 
of changing the time interval between successive 
accretion events. We have tested, for accretion 
amounts of O.OSMjup, the cases of events separated 
by 1 Myr and 4 Myr. The results are presented 
in Figure [12] for comparison with Figure [S] Figure 
[13] shows the surface lithium abundance variation 
with time, together with the previously studied 
case of 2 Myr. The left graph displays the results 
for five events and the right graph for ten events. 
In each case, two series of models are computed : 
one with Ct=12 and the other one with Ct=1000. 
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Fig. 7. — Lithium surface abundances over the accretion/mixing period in models experiencing 5 accretion 
episodes of O.OSMjup. The presented models have different initial metallicities ([Fc/H]=-0.20, 0.00, 0.16, 
0.20, 0.24) and different masses (KRT coefficient, Cf=:12). For the I.2OM0 models, the five curves merge. 
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The final results are not trivial. For larger time 
intervals, the /.t-gradients before every new accre- 
tion event are smaller, so that the diffusion coef- 
ficients and thus the mixing efficiencies are also 
smaller, leading to less lithium destruction. On 
the other hand, the total mixing episodes last 
longer, in favor of more lithium destruction. Fig- 
ure [13] shows that these two combined effects can 
lead to opposite results according to the circum- 
stances. For Ct=12, the results obtained with the 
three different time intervals are very similar, al- 
though lithium is slightly less depleted for increas- 
ing intervals. For Ct=1000, the situation is differ- 
ent if accretion stops after five or ten events. For 
five events, the results are also very similar but this 
time the case of events separated by 1 Myr lead to 
slightly less depletion than the other cases. On the 
other hand, for ten events, the situation is similar 
to the case described for Ct — 12, except that the 
difference is more important. An interesting point 
in this case is that the end of the accretion pe- 
riod is clearly visible in the graph. It corresponds 
to the time when no fresh lithium is brought any- 
more, so that the lithium abundance variation is 
only due to nuclear destruction. 
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Fig. 9. — Same figure as Figure [3] for a I.IOM© models experiencing 10 accretion events of O.OSMjup, 
computed with Ct=1000. The Hue convention is the same as on Figure [3] except for the hght magenta 
dashed hues which represent the profiles at 30 and 51 Myr on the ZAMS (from top to bottom). The black 
arrow on the upper right panel shows the extension of the mixed zone after 10 accretion episodes. A small 
stable helium zone develops after one accretion event and never come back afterwards. 
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in Figure O 
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Fig. 10. — Same figure as Figurc|4]for models com- 
puted with Ct = 1000. Dashed lines are models 
experiencing (from bottom to top) : 1, 5 or 10 
accretion events of IMg. Solid lines are models 
experiencing (from bottom to top) : 1, 5 or 10 
accretion events of O.OSMjup. 



5.7. Computations using recent diffusion 
coefficients deduced from numerical 
simulations 

As already mentionned, the expressions for the 
diffusion coefficients deduced from recent numer- 
ical simulations lead to results close to the KRT 
ones with Ct = 12, except for the boundaries of 
the thermohaline regime. Here we are basically 
interested in the Rq = 1 limit, which corresponds 
to the boundary between the dynamical convec- 
tion and thermohaline mixing. In the KRT pre- 
scription, for which Dth becomes very large at this 
limit, there is a simple continuity between the two 
regimes. After accretion, if i?o < I7 the region be- 
low the standard convectivc zone becomes dynami- 
cally unstable, which leads to complete mixing and 
may be modelled by a very large diffusion coeffi- 
cient. As the ^-gradient decreases, i?o increases. 
When it reaches unity, thermohaline convection 
begins. One thus expects a continuity between 
the diffusion coefficients, which means that when 
Rq becomes slightly larger than one, Dth should 
rapidly decrease from a very large value to the 
value obtained in the intermediate regime. This 
description is well fitted by the KRT coefficient, 
but not by the TGS one where Dth remains finite 
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Fig. 12. — Molecular weight and diffusion coefficients during the accretion/mixing process in two I.IOMq 
models experiencing 10 accretion episodes of O.OSMjup with timcstcp of 1 or 4 Myr. The black lines represent 
the /x-profiles at 2Myr on the ZAMS, just before the first accretion event. The dark blue to magenta solid 
lines present the profiles after each new accretion event. The light magenta dashed lines show profiles when 
the accretion has stopped but when the thermohaline mixing still proceeds : at 26, 107 and 229 Myr for the 
left panel, at 57, 122 and 260 Myr for the right panel. The horizontal arrows in the lower panels locate the 
mixed region after 10 accretion events. 
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Fig. 13. — Lithium depletion in models computed with different time intervals {At) between successive 
accretion events. Left figure : models experiencing 5 accretion events of 0.03Mjup. Right figure : models 
experiencing 10 accretion events of O.OSMjup. For each accretion scenario, two series of models are computed 
: one with Ct=12, one with Ct=1000. For each Ct value three models are presented assuming different At 
values : 1 Myr (dashed-dotted lines), 2 Myr (solid lines), 4 Myr (dashed lines). 
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as it reaches unity. However, a close look at their 
Figure 2, which presents the turbulent heat and 
compositional fluxes as a function of their rescaled 
density ratio r, reveals a clear tendency for the nu- 
merical simulation points to go up to large values 
for r — 0, corresponding to i?o = 1- We thus infer 
that their should be a small transition region be- 
tween the TGS value and the very large one needed 
for i?o < 1, consistently with the simulations. On 
the other hand, the value of the mixing coefficient 
as given by TGS drops to zero at the other limit, 
for r = 1, corresponding to i?o = I/t". This fea- 
ture, which is not present in the KRT coefficient, 
is consistent with the fact that Dth must become 
zero for Rq > 1/r. 

Figure [H displays the KRT and the TGS co- 
efficients precisely computed as a function of Rq 
in one of our models (M=1.1OM0, [Fe/H]==0.20 
at 4 Myr). T^^^ gr aph is similar to Figure 3 of 



Traxler et al.l (|201l[ ). The difference comes from 



the fact that TGS used constant k^, kt and 
in their computations, whereas we use the exact 
values computed in our model. These values vary 
with depth and with age, as well as the /i-gradient. 

The two coefRcients are very similar in the re- 
gion between i?o = 5 x 10^ and Rq ~ 5 x 10^, 
i.e. between l/50r and 1/5t, with r = 2.5 x 10^. 
For Rq close to one, Dtgs has a finite value 
which would produce a discontinuity with the case 
i?o < 1. For this reason, we have added a transi- 
tion region in our computations (solid black line, 
label Dt). For the limit where Rq becomes close 
to 1/r, the TGS coefficient decreases rapidly, con- 
trary to the KRT one. Also shown on this graph i s 
the coefficient used bv lCharbonnel fc ZahnI ( 20071 ). 



labelled Dcz- 

Here we present results obtained with the TGS 
diffusion coefficient, with a very small transition 
region below r = 0.02, similar to the KRT one, to 
preserve the continuity with the dynamical regime. 
Figure[T5]displays the lithium and beryllium abun- 
dance variations with time for various accretion 
scenarios, which must be compared with Figure 5. 
We can see that the obtained lithium and beryl- 
lium depletions are very similar, although slightly 
smaller, that those obtained with the KRT pre- 
scription with Ci = 12. In the worst case, the 
relative difference is of order 25%. On the con- 
trary, the depletions obtained with Ct = 1000 are 
much too large. 
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Fig. 14. — Thermohaline diffusion coefficients 
(cm^.s~^) computed in a I.IOM0 model with 
[Fe/H]=0.20 at 4 Myr. Long dashed-small dashed 
(green) line: KRT coefficient; dashed (blue) 
line: TGS coefficient; sofid (black) fine: TGS 
coefficient including a transition region for Rq 
close to one; dotted line: coefficient used by 



Charbonnel fc ZahnI (120071 ). 
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Fig. 15. — Same as Figures 5 with the lTraxler at al.l (|2011l ) prescription for the mixing coefficient. 



6. Summary and discussion 

We have investigated the effects of planetesimal 
accretion onto solar-type stars. We have shown 
that the dilution of metal-rich material in the 
surface convective zone of young main sequence 
stars creates an unstable //-gradient at the tran- 
sition between the convective and radiative zones. 
The induced double diffusive instabilities, the so- 
called thermohaline convection, result in a mixing 
which dilutes the metal-rich matter until the fi- 
stratification becomes stable. The thermohaline 
convection effects combined with the stabilizing 
action of atomic diffusion softens the inverse /i- 
gradient on a few millions or tens of million years, 
depending on the considered accretion scenario. 
Note that this ti me scale i s much larger than the 
one suggested by IVauclaiil (|2004h . The reason is 
that, following KRT, she used as a characteristic 
length the evaluated fingers dimension, whereas 
the real length to be used is that of the whole 
mixed region. At the end of the accretion/mixing 
period, only a very small metallicity increase may 
remain, much smaller than the average overmetal- 
licity observed in exoplanet-host stars. This re- 
sult is consistent with a primordial origin of this 
overmetallicity. On the other hand, as suggested 



by iGaraud ( 2011 ). the thermohaline mixing al- 
lows connecting the stellar surface with the Li- 
nuclear burning region, which may lead to rapid 
Li-depletion depending on the accretion rate and 
the number of impacts experienced by the star. 

The Li destruction depends on many different 
parameters linked to the stellar structure and to 
the rate, frequency and amount of the accretion 
episodes. Concerning the stellar structure, the 
lithium depletion is obviously larger for deeper 
convective zones, which is related to the stellar 
mass and metallicity. It also depends on the age at 
which the accretion occurs, because the diffusion- 
induced /i-gradients act as stabilizers for the ther- 
mohaline convection. Concerning the accretion 
process itself, wc have shown that different re- 
sults are obtained not only for different accreted 
amounts, but also according to the way it is ac- 
creted : all at ones, or in successive episodes. 
Different stellar characteristics and various accre- 
tion scenarios produce, at the end of the accre- 
tion/mixing period, a significant dispersion in the 
Li abundances. 

The results also depend on the computational 
parameters used for the treatment of convection. 
In the present computations, the convective zone 
is computed whithin the framework of the mixing 
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length theory, with a mixing length ratio adjusted 
on the solar models {a = 1.8). No overshooting 
nor extra-mixing except thermohalinc convection 
is added. Thermohaline convection is treated us- 
ing the diffusion approximation. 

We have tested several prescriptions for the dif- 
fusion coefficient. Recent 3D numerical simula- 



tions bv lTraxler et al.l (j201lh allowed real improve- 
ments in the modelling. We have shown that the 
corresponding prescription for the diffusion coef- 
ficient leads to results very similar to the KRT 
one, provided that the aspect ratio to the fingers, 
that is the ratio of their length to their width, be 
taken equal to one as suggested by the simulations 
(see also ,Dcnisscnkov (2010^). This corrcsonds to 
a numerical coefficient Ct = 12 in the KRT ex- 
pression. As the diffusion coefficients arc larger 
for larger /i gradients, which themselves decrease 
faster when the diffusion processes are stronger, 
the final lithium depletion is not simply related 
to the initial diffusion coefficient value. Here we 
have tested three cases: the KRT coefficient with 
Ct = 12, the same coefficient with Ct = 1000, as 
suggested bv ICharbonnel &: ZahnI ( 2007 ) for the 
RGBs and the coefficient deduced from numerical 
simulations by TGS. We confirmed that the KRT, 
Ct = 12 coefficient gives results similar to the TGS 
one. 

Garaudl(l201lh performed computations of ther- 
mohaline convection in a IAMq static model, us- 
ing the TGS coefficient. No atomic diffusion was 
introduced in her computations. We can see from 
her results, and particularly from the values of 
the diffusion coefficient and the mean molecular 
weigth presented in her Figure 3, that our results 
and hers are similar, of the same order of magni- 
tude. This convergence is quite encouraging. 

The general trend of lithium depletion with age 
in various galactic clusters may be accounted for 
by rotation- induced mixing and/or internal waves 



in a som ewhat satisfactory way fe.glTheado fc Vauclaii 
(l2003bl ): iTalon fc Charbo^ ^00^). The large 



dispersion observed in stars of similar intrinsic pa- 
rameters like mass, age, rotation velocity, chemi- 
cal composition nee ds a more complicate process 



to be accounted for (ISqusa et al.ll2010l: iKing et al 



2010l ) . As discussed by Pinsonneault et al. ( 2011 ) , 
the details of the rotation history of stars can help 
explaining these different features. Here we sug- 
gest that the accretion/thermohaline mixing pro- 



cess may also be an important parameter in the 
description of the lithium dispersion in solar-type 
stars, and that it should be taken into account. 

Comparing our results with the observed abun- 
dances of lithium and beryllium in solar type stars 
lead to several concluding point. 

1) If one assumes that accretion of planetary 
matter onto solar type stars is a common phe- 
nomenon, thermohalinc convection should often 
occur. On the other hand, if most of the lithium 
abundance behavior in solar type stars is well ac- 
counted for by rotation-induced mixing and re- 
lated processes, as generally believed, it leads to 
constraints on the accretion efficiency. Using the 
realistic TGS or KRT thermohaline coefficients, 
we find that the stars must not accrete more than 
ten earth masses after their arrival on the Main Se- 
quence, to keep their original lithium value within 
0.1 dex. 

2) Stars which accrete more than ten earth 
masses are expected to be subject to stronger 
lithium deplet i on. A s alr eady po i nted o ut by 



Theado et al. 



(|20i3) and iGaraudI (|201lh . this 
process could be an elegant explanation for the ob- 
servations of stars in which strong lithium deple- 
tion is found. Note that this does not only apply to 
those stars for which planets have been detected, 
as stars may well host planets even if they have not 
yet been detected because of observational bias. 

3) In any case, the observed ovcrmetallicity in 
exoplanets-host stars is not due to accretion, and 
it is not possible either to explain any increase of 
lithium or of the ^Li/^Li ratio by accretion, be- 
cause of the thermohalinc dilution. 

4) The solar case deserves a more precise anal- 
ysis. Whether accretion at the period of the late 
bombardment has participated in the presently 
observed lithium depletion or not depends on the 
amount of matter swallowed by the Sun, and also 
on the diffusion-induced helium gradient at the 
time of the event. 

5) Finally, a complete study of lithium deple- 
tion in solar-type stars needs a combined computa- 
tion of this process with other depiction processes 
like rotation-induced mixing. Complete evolution- 
ary sequences including accretion, thermohaline 
mixing and rotational mixing will be investigated 
and presented in a forthcoming paper. 
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